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Abstract 

In the present study, we establish the existence of nontrivial site percolation 
threshold in the Relative Neighborhood Graph (RNG) for Poisson stationary point 
process with unit intensity in the plane. 

1 Introduction 

Percolation theory is very useful to describe various physical phenomena. In particular, 
there are important connections with phase transition problems [15, 14J. 

The interest for percolation problems has grown rapidly during the last decades: see 
Lyons and Peres fT9l for percolation on trees and networks and Meester and Roy ll20l 
for continuum percolation and the references therein. In 1996, Haggstrom and Meester 
IfTTl proposed results for continuum percolation problems for the fc-nearest neighbor 
graph under Poisson process. In a recent paper, Balister and Bollobas f3l give bounds 
on k for the fc-nearest neighbor graph for percolation with several possible definitions. 

Site and bond percolation as well as phase transition for several statistical me- 
chanics models have also been studied by Haggstrom [T6| for general graphs and in 
particular for the Delaunay triangulation for the Poisson point process. It also shows 
for graphs with a bounded degree if the bond percolation threshold is not trivial then 
it implies a phase transition for the Ising model on this graph. For general results on 
Delaunay graphs and Voronoi tesselations see Moller 1,21] . For recent results on perco- 
lation on these graphs see Balister et al. ||4]|2l. 

Benjamini and Schramm ||6l proposed a comprehensive study on general graphs, 
with special focus on Cayley graphs, quasi-transitive graphs and planar graphs. Re- 
cently, Procacci and Scoppola f22] proposed sufficient conditions on infinite graphs 
to deduce a non trivial bond percolation threshold. Among these assumptions, they 
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assume that the dual graph is bounded degree. It is interesting to relax this condi- 
tion in order to deal with proximity random graphs (which in general have a dual not 
bounded degree) like the skeletons on point Poisson processes of the plane. These ran- 
dom graphs includes the Gabriel graph and the Relative Neighborhood Graph (RNG) 
which are important for many applications. 

Remark that continuous models defined on nearest neighbors graphs are interesting 
for small temperature as alternative of standard models on regular networks, because 
it allows vibrations and deformations of the network and may be find an application in 
physics of the solid state. In particular one example is the study of the order-disorder 
transition of binary alloys or ionic cristals. It is well-known that Delaunay graph or 
Voronoi regions (rather called Wigner-Seitz grid and Brillouin zone in physics frame- 
work) take a fondamental place for the understanding of the electrical current, waves 
propagation and phase transitions observed by Bragg diffraction of X rays. 

Another domain of application should be found in cancerology for the study of 
the growth of tumour when the cancer cells suddently begin to invade healthy tissue. 
The Delaunay graph is well adapted for such study as explained in |10| but the RNG 
should give some more information. More precisely, in histology we have some slides 
with marked cells (cancer cells and normal cells in first approximation) and probably, a 
connection between percolation and the grade of a given cancer will be helpful to give 
an aid for the diagnostic of a pathologist. 

The existence and unicity at small activity of nearest neighbors stationary Gibbs 
states can be found in Bertin et al. |7 8, 9|. In Bertin et al. lfT2l . the phase transition 
in the Delaunay continuum Potts model is established. It is a generalization of the 
Lebowitz-Lieb model as described in Georgii and Haggstrom |13| where the soft re- 
pulsion between several species of particules acts on the Delaunay graph. What is the 
good Delaunay subgraph on which the repulsion is strong enough to maintain a phase 
transition? In terms of percolation, it means: is bond percolation maintained in this 
subgraph? Bertin et al. Ill] gave an answer for the Gabriel graph. Another well-known 
subgraph of Delaunay graph is the Minimum Spanning Tree (MST). The structure of 
spanning forest generalized minimum spanning tree for infinite graphs. Constructed 
with the greedy algorithm, this graph exactly is described in the book of Meester and 
Roy II20I . The link between branching number and percolation on trees is proved by 
Lyons ifTSl . Connection between minimum spanning forest and occupied and vacant 
percolation is strong (for more details on simultaneous uniqueness see Alexander |T|). 
Thus, the minimum spanning forest is a tree with one infinite path a.s. in two dimen- 
sions when the points are distributed under a stationary Poisson process. It comes that 
a.s. the site or bond critical threshold are equal to 1. 

The present study gives an answer for the RNG, well known in computational ge- 
ometry, when the points are distributed under a stationary Poisson point process with 
unit intensity in the plane. We adapted a powerful method of the rolling ball proposed 
by Balister and Bollobas [SJ relying on 1 -independent bond percolation on 1? . Then 
if we control the probability of having less than a fixed number of points in a given re- 
gion and considering the event that all the sites are open in this region, we can proceed 
similarly as in Haggstrom l,16il . 

The paper is organized as follows. The first section is devoted to some definitions 
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and notations. Next, the main result on site and bond percolation on the RNG is given. 
Then, we give the proof of the main result, by using a result of bond percolation in 
the 1 -independent case in 1? and the rolling ball method. We conclude on possible 
extensions of this work. 



2 Notations and definitions 

Let \A\ denotes the Lebesgue measure if the set A is a bounded Borel set of R^, and 
the counting measure if A is a discrete set. 

Given a finite box A C K^, we denote by Ha the Poisson point process on the 
locally finite set of points in A denoted by JIa with intensity 1 i.e. 

/ /dllA = exp(-|A|) ^ / /({xi, . . .x„})(ia;i, . . .(ia;„ 

for any bounded measurable function f on JIa- 

Let the set of locally finite subsets of K^. We have to consider only the config- 
urations e which are in general position (four points on the same circle do not 
occur and no three points are colinear) in order to ensure the existence and unicity of 
the Delaunay graph. One can notice that, for any stationary point processes, the proba- 
bility of the set of tesselations in general position is equal to one 1.2 1 J . Let us recall the 
definitions of Delaunay, Gabriel and relative neighborhood graphs. 

Definition 1 The Delaunay graph Deli (}p) of a configuration ip in f2 is the set of edges 
of the unique triangulation Del^{Lp) in which the interior of the circle circumscribed 
by every triangle of Del^{ip) does not contain any point ofip. 

Definition 2 The Gabriel graph Gab{ip) of a configuration ip in Q is defined as the set 
of edges {u, v} d ip such that the open circle with {u, v} as diameter does not contain 
any point of the configuration ip. 

Definitions The Relative Neighborhood graph Rn(ip) of a configuration (p in f7 is 
defined as the set of edges {u, v} <Z such that the intersection of the disks with 
center u and v with radius uv does not contain any point of the configuration (p. 

These graphs are planar in R^. Furthermore, the RNG is a subgraph of the Gabriel 
graph which is a subgraph of the Delaunay Graph (see figure[T]i. 



3 Main result 

We first deal with the site percolation on the RNG for the Poisson point process 
We introduce the Bernoulli process 8(<i>) providing the type picking mechanism (1 for 
open and for closed) of the points (or sites) in Obviously, 4> = ($, 8(<i>)) can be 
seen as a marked Poisson process. The probability measure of $ is given by : 
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Figure 1 : Vacuity regions for Gabriel and Relative neighborhood graphs for an edge 
uv 

where /i^ is the probability measure on {0, 1}"^ of the Bernoulli process Q{(p) given 

the configuration ip ^ ft. Similarly is the probability measure on 

of the Bernoulli process @{Rn{(p)) given the graph Rn{(p) where ip E fl. We define 

pf*%ip, Rn{ip)) andp^°"''((^,i?n(v?)) such that: 

/i^(3 at least one infinite open cluster in Rn{ip)) 

_ r 1 if p> pf^%f,Rn{ip)) 
^ \ if pKpf iip^Rniip)) 

l^'Rn{ip) '^^ least one infinite open cluster in Rn{ip)) 

_ f 1 if p> p''°'"^{ip,Rn{Lp)) 
~ \ if p <pl°'^'^{ip,Rn{ip)) 

Let us recall fTSl the following well known relation between p^***^ and p^""*^ on a given 
graph with bounded degree: as the graph [tp^ Rn{ip)) have a degree at most 6, 

y^en, l/5<p''r''{^,Rn{p))<pf'^{ip,Rn{p)) < 1 - [1 V ^"M)]'- 

We now introduce {Rn, 11) defined as the lowest p for which the probability 
of the event that there exists an infinite open cluster in the RNG relative to the marked 
point Poisson process $ is equal to 1. By ergodicity of 11, this previous event, invariant 
by translation, is a trivial event. The marked point Poisson process ^> may exhibit 
some percolation phenomenon with critical value j<^''° {Rn, 11). We want to prove in 
the following theorem that this is a non trivial critical value, i.e.. 

Theorem 1 {Rn, U) < 1 andp^°'"^ {Rn, U) < 1. 

The following section is devoted to the proof of this theorem. 
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4 Proof 



First, we point out that the method proposed by Haggstrom lfT6l and used for k nearest 
neighbor graph and Delaunay graph, as well as the adaptation by Bertin et al. ifTTIl to 
the Gabriel graph does not applied to the RNG. However, it is not sufficient to control 
the probability of having at least one point and less than a fixed number of points in 
each small box K (see figure|2]l. 

Indeed, we may choose some configuration of points such that the length and num- 
ber of points of a Rn— path between two points are arbitrary large, see figure |2l wt 
is an edge of the RNG but uv is not because for example the point q belongs to the 
vacuity region of this edge. 



Figure 2: Example of construction of configuration with arbitrary large Rn- path. 

Moreover, such methods are based on comparison with independent bond percola- 
tion on the grid Z^. As, in Balister and Bollobas Q in the case of k nearest neighbor 
graph, we adapt the method of rolling ball (see figure |3]l to the RNG. We procceed in 
two steps: 

1 . Controhng the probabihty of some suitable configurations of points under Pois- 
son point process. 

2. Consider the Bernoulli site percolation in such configurations. 

To prove that continuous percolation occurs, we shall compare the process to var- 
ious bond percolation models on I?. In these models, the states of the edges will not 
be independent. However they will satisfy the following definition: 

Definition 4 A bond percolation model is 1 -independent if whenever Ei and E2 are 
sets of edges at graph distance at least 1 from each another (i.e., if no edge of Ei is 
incident to an edge of E2) then the state of the edges in Ei is independent of the state 
of the edges in E2. 

We shall use the following result in Balister et al. ||5]. 
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Figure 3: The Rolling Ball Method 



Theorem 2 If every edge in a 1 -independent bond percolation model on 1? is open 
with probability at least 0.8639, then almost surely there is an infinite open compo- 
nent. Moreover, for any bounded region, there is almost surely a cycle of open edges 
surrounding this region. 

Let us first consider the case of percolation in the RNG. Write u ^ v \f uv \s an 
edge of the underlying graph Rn{(p). For percolation we need to find an infinite path, 
i.e., a sequence ui,U2 ■ . ■ with Ui ~ u^+i for all i. Consider the rectangular region 
consisting of two adjacent squares Si, S2 shown in figure |3] Both Si and ^2 have side 
length 2r + 2s, where r and s are to be chosen later We define the basic good event 
£si,S2 to be the event that every vertex ui in the central disk Ci of Si is joined to at 
least one vertex v in the central disk C2 of ^2 by a Rn— path, regardless of the state of 
the Poisson process outside of Si and ^2. 

Now consider the following percolation model on Z^. Each vertex E 1? 

corresponds to a square \Ri, R(i + 1)] x [_Rj, R{i + 1)] e M^, where R ~ 2r + 2s, 
and an edge is open between adjacent vertices (corresponding to squares 5*1 and S2) 
if both the corresponding basic good events £si.S2 ™d £32, Si hold. Note that this 
is indeed a 1 -independent model on Z'^ since the event £si,S2 depends only on the 
Poisson process within the region Si and 6*2, and thus sets of edges at distance at least 
one apart in 1? depend on the Poisson process in disjoint regions of R^. Any open 
path in 1? corresponds to a sequence of basic good events £^5^ 5^ , £32,83 ■ ■ ■ that occur, 
where Si is the square associated with a site in 7?. Every vertex ui of the original 
Poisson process that lies in the central disk Ci of Si now has an infinite path leading 
away from it, since one can find points Ui in the central disk of 5*^ and paths from 
Ui-i to Ui inductively for every i > 1. In particular, each such ui lies in an infinite 
component. Moreover, such vertices exist in Ci, so there is an infinite component. One 
can choose r and s so that the probability that the intersection of a basic good event is 
large and then we will apply the theorem |2l 
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In order to bound the probability that this intersection of a basic good event fails, 
we shall use the following rolling ball method. Let Ci, C2, and L be as in Figure[3] {L 
is the region between the two disks Ci and C2.) We need to define Es^ ^^ the event 
that for every point u G Ci U L, there is a u such that: 

a) ti w; 

b) d{u, v) < s; and 

c) u G Dy, where D„ is the disk of radius r inside Ci U L U C2 with v on its Ci- 
side boundary (the dotted disk in Figure[3]l. Note in particular that (b) implies that the 
condition u ^ v in (a) is independent of the Poisson process outside of Si U 82- This 
is because both u and v are at distance at least s from the exterior of Si U 5*2, so the 
event that {u, v) is an edge of the RNG only depends on the points within U We 
denote Eis-^.S2 the complementary of E5j^,s2- The probability of £5^ 53 is bounded by 
the expected number of points u for which above conditions (a)-(c) fail. Thus, we have 

n(E5i,5j < 2r(2r + 2s)pRn.r,s (1) 

where pRn.r,s is the probability that (a)-(c) fail for some fixed v. Notice that this 
probability is independent of the location of w in Ci U L. 

Lemma 1 We can choose r and s such that T\{EiSi,S2) arbitrary small. 

Proof: 

Let a) be the disk of radius a and of center x. 



Then, 

Pfl„,r,. <e-|^"^^('''^)l+ / Pr.s{u)du 

JD^nD{v,s) 

wherepr,8{u) = e-l^"nD(«, d(t.,,;))| _ ^-\D{v,d{u,v))nD{u,diu,v})\D^.\^ the prob- 
ability that u is the closest point to v inside Dy, but that {u, v) is not an edge of the 
RNG. To calculate this upper bound, note that 

\Dy n D{v,s)\ = -rs^ 1 - ■^+(2r^ ~ s2)arcsin (^-^^ + sV/2. 
By choosing polar coordinates (a, 6) of u, it comes: 



' D„nD{v,s) 



arccos(Q/(2r)) 



J(a, e) de da 



where 

J{a,e) = 1 - e-\D{v,a)nD({afi),a)\D^\^ 

To calculate this last integration, we have : 

^arccos(^)-7r/3 

J(a, 9) d9 = arccos — 7r/3 - 



/•arccos^ 
^ arccos^ ^ ^ —7r/3 



J{a,9)d9 = 7r/3- 2- 



^—\L{a,r, arccos ^^^—77/3) I 



a-" 



1 - e 



1 



D{v, a) 



D(u, a) 



Figure 4: The Lune L{a, r, 9) = D(u, a)\Dy is the dotted area 



with 



|L(a, r, 9)\ = c?d + {a^ — r^) arcsin 



asmf 



. V?'^ + — 2ar cos d 
|L(a,r,arccos(f ) -7r/3)| = ^(f - ^) -r^ arcsin ( 



+ ar sin i 



\1r) 



4^2 • 



Thus, 



/ Vr s{u) du =2 r ae-l^''"^('''")l x (arccos f ^) 

JD^nD{v,s) ' Jo I \2r/ 

f 

Jo 



g-|Z,(a,r,e)| _ _^ 



Take 



/ Pr,s{u)du =2 Z'aarccosf;^) e-l^-^^f^'^^lrfa 

JD^nD(v,s) ' Jo \2rJ 



-2 rae-l^"nD(.,a)| /■ 
/o„ ^0 



arccos (^)—7r/3 



_4 / lg-|D„nD(i;,a)|-|I,(a,r,arccos(^)-7r/3)| 

a 



da. 



e-\i^M)\deda 



da 



we conclude that 

PRn,r,s < e-|^"nD(.,.)| Z 

>/£), 



1 



Pr,s(u) du 



Jo OL 



vn£)(t>,s) 

/■s /■a 



-|£>„nD(ti,a)|-|I,(a,r,arccos(^)-7r/3)| 



da 
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which can be bounded by 



PRn, r,s < 1-2 / ae 



-|-D„nD(i>,a)|-|L(Q,r,arccos(^)-7r/3)| 



arccos (^) - tt/S + ^^^—^ 



da. 

□ 



For instance, the bound involved in inequahty[T]gives around 10 '^'^ with r — s = 
8000. 

For our purpose, we also need to control the probability for the Poisson point pro- 
cess of having at least one point in Ci and less than m points in C2 U Ci U L. 

We denote Fc, = {*(Ci) > 1} and = {$(6*2 U Ci U L) < m}. Notice that 
as Es,.S2 n Fc, n Am C £si,S2 then 

Esi,s, n Es,,5, n Fc, n Fc, n a„ c £si,s2 n £s2.s^■ 

We obtain that 

n(fsi.s. nfs.A) > 1- [n(Es,,sJ + n(E5,.sJ + n(FcJ + n(FcJ + n(A„o] 

where n(Fci) = UiFc^) = e"'"'' and 

/i;>m 

Choosing r — s, we have the following bound 

n(A„)< "7''T' ^ 

(m + 1)! 

This bound becomes negligible whenever r ^ 1 and m > e(8 + 7r)r^ using Stirling 
formula. 

Set e = 0.1361. Similarly as in Haggstrom lfT6l but adapted in the 1-independent 
case, let Br^s be the event that all the sites are open in C2 U Ci U L with probability 
p = I — ^ and also define Cr,s — £si,S2 ^ ^82,81 H S^.s- Then, we have: 



JEsj nEs2 ,si nJ^ci nJ^ca nA„ 

> (1 - e/2y" > 1 - e ^ 0.8639 
because we can choose r, s, m (as preceding) such that 

n(Es,,s2 n Es,,si n Fc, n ^^c^ n A™) > i - e/2 

We conclude with theorem|2]that W) < I - ^ < I. 
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5 Concluding remarks 



This kind of proof also apply in the case of bond or site percolation of the k nearest 
neighbor graph: it is direct consequence of |3 |. It is sufficient to use what they called 
Pu, PB, Pi or po for several possible definitions of percolation. 

As suggested in |[3l with a high confidence, fc = 3 is the critical out-degrees for 
percolation on the k nearest neighbor graph. We notice that in the RNG the number 
of neighbors is bounded by 6. So a point have neighbors in several directions. It may 
be interesting to study a family of Delaunay subgraphs defined on sequence of vacuity 
regions such that the number of neighbors of each point is the lowest as possible but 
keep good connectivity properties for percolation purposes. Otherwise a challenge 
would be to extend the method of the rolling ball when points are distributed under a 
Gibbs point process for instance a hard-core point process. 
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